The chiral transitions in heavy-light mesons 

A.M.BadalianO Yu.A.SimonovQ and M.A.Truso\! 
ITEP, Moscow, Russia 
(Dated: December 22, 2007) 

Abstract 

The mass shifts of the P-wave D s and B s mesons due to coupling to DK, D*K and BK, 
B*K channels are studied using the chiral quark-pion Lagrangian without fitting parameters. The 
strong mass shifts down ~ 140 MeV and ~ 100 MeV for D*(0 + ) and D s (l + ') and ~ 100 MeV for 
B*(0+) and B s (l+ ) are calculated. Two factors are essential for large mass shifts: strong coupling 
of the + and 1 + ' states to the S-wave decay channel, containing a Nambu-Goldstone meson, 
and the chiral flip transitions due to the bispinor structure of both heavy- light mesons. The 
masses M(B*(0+)) = 5695(10) MeV and M(B S (1+')) = 5730(15) MeV,very close to Af(B(0+)) 
and M(B(1 + ')), are predicted. Experimental limit on the width r(D s i(2536)) < 2.3 MeV puts 
strong restrictions on admittable mixing angle between the 1 + and 1 + ' states, \(f>\ < 6°, which 
corresponds to the mixing ang le 6 between the 3 Pi and 1 P 1 states, 29° < 6 < 41°. 
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I. INTRODUCTION 



The heavy-light (HL) mesons play a special role in hadron spectroscopy. First of all, a 
HL meson is the simplest system, containing one light quark in the field of almost static 
heavy antiquark, and that allows to study quark (meson) chiral properties. The discovery 
of the .D s (2317) and -D s (2460) mesons with surprisingly small widths and low masses 

has given an important impetus to study chiral dynamics and raised the question why their 
masses are considerably lower than expected values in different approaches: in relativistic 
quark model calculations |^|-B|, on the lattice p|, in QCD Sum Rules (si, [sj], in chiral models 
10 ] — 12 1 ( for reviews see also [la. \h\ ) . The masses of D s (0 + ) and D s (l + ') in closed-channel 



approximation typically exceed by ~ 140 and 90 MeV their experimental numbers. 

Thus main theoretical goal is to understand dynamical mechanism responsible for such 
large mass shifts of the + and 1 + ' levels (both states have the light quark orbital angular 
momentum I = and j = 1/2) and explain why the position of other two levels (with 
j = 3/2) remains practically unchanged. The importance of second fact has been underlined 



by S.Godfrey in 



5]. 



The mass shifts of the -D s (0 + , 1 + ) mesons have already been considered in a number 
of papers with the use of unitarized coupled-channel model 151, in nonrelativistic Cornell 

n nn 

model [16|], and in different chiral models [17|-[19[]. Here we address again this problem with 
the aim to calculate also the mass shifts of the D s (l + ') and B S (Q+, 1 + ') states and the widths 
of the 2 + and 1 + states, following the approach developed in 18J, for which strong coupling 
to the S-wave decay channel, containing a pseudoscalar (P) Nambu-Goldstone (NG) meson, 
is crucially important. Therefore in this approach principal difference exists between vector- 
vector iVV) and VP (or PP) channels. This analysis of two-channel system is performed 
with the use of the chiral quark-pion Lagrangian which has been derived directly from the 
QCD Lagrangian [^(J and does not contain fitting parameters, so that the shift of the D*(0 + ) 
state ~ 140 MeV is only determined by the conventional decay constant fx- 

Here the term "chiral dynamics" implies the mechanism by which in the transition from 
one HL meson to another the octet of the NG mesons <fi is emitted. The corresponding 
Lagrangian AL FGM , 

AL FCM = q(er) exp(z7 5 0/ f w )q, (1) 
contains the light-quark part exp(z7 5 0//^), where is the SU(3) octet of NG mesons and 
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the important factor 75 is present. In the lowest order in ei> this Lagrangian coincides with 
well-known effective Lagrangian AL e ff suggested in 21], [22], where, however, an arbitrary 



constant g& is introduced . At large N c , sued in 



unity, g A = 1. In 



[lO, 



21 
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this constant has to be equal 
22] this effective Lagrangian was applied to describe decays of HL 



mesons taking g& < 0.80. 

More general Lag rangian AL fcm CQ) was derived in the framework of the field correlator 



method (FCM) 20j, |23j , in which the constant g A = 1 in all cases, and which contains NG 



mesons to all orders, as seen from its explicit expression (TjQ). 

In Appendix A with the use of the Dirac equation we show that in the lowest order in 
cf) AL F cm = AL e ff, if indeed g A — 1. In our calculations we always use AL FC m with the 
gA — 1 and derive the nonlinear equation for the energy shift and width, AE = AE — ^, 
as in 18|. We do not assume any chiral dynamics for the unperturbed levels, which are 
calculated here with the use of the QCD string Hamiltonian 2J, |25[ , because the mass shift 
AE appears to be weakly dependent on the position of unperturbed level. 

It is essential that resulting shifts of the J p (0 + , 1 + ') levels are large only for the D S ,B S 
mesons, which lie close to the DK, D*K, BK, B*K thresholds, but not for the D(1P), B(IP) 
mesons, in this way violating symmetry between them (this symmetry is possible in close- 
channel approximation). In our calculations shifted masses of the D s (0 + ) and B s (0 + ) practi- 
cally coincide with those for the D(0 + ) and B(0 + ), in agreement with the experimental fact 
that M cxp ( J D(0 + )) = 2350 ± 50 MeV ^ is equal or even larger than M exp (D s (0 + )) = 2317 
MeV. The states with j = 3/2 D s (l + , 2 + ) and -D(l + , 2 + ) have no mass shifts and for them 
the mass difference is ~ 100 MeV, that just corresponds to the mass difference between the 
s and light quark dynamical masses. 

For the D s (l + ') and B s (l + ') mesons calculated masses are also close to those of the D and 
B mesons. Therefore for given chiral dynamics the J p (0 + , 1 + ) states cannot be considered 



lit 



as the chiral partners of the ground-state multiplet J p (0~, 1~), as suggested in 

We also analyse why two other members of the IP multiplet, with J p = 2 + and 1 + , do 
not acquire the mass shifts due to decay channel coupling (DCC) and have small widths. 
Such situation occurs if the states 1 + and 1 + ' appear to be almost pure J = § and j — \ 



states. Still small mixing angle between them, |0| < 6°, is shown to be compatible with 
experimental restriction on the width of -D s i(2536), admitting possible admixture of other 
component in the wave function (w.f.) < 10%. 
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In our analysis the 4-component (Dirac) structure of the light quark w.f. is crucially 
important. Specifically, the emission of a NG meson is accompanied with the 75 factor 
which permutes higher and lower components of the Dirac bispinors. For the j = 1/2, P 
-wave and the j = 1/2, S -wave states it is exactly the case that this "permuted overlap" of 
the w.f. is maximal because the lower component of the first state is similar to the higher 
component of the second state and vice versa. We do not know other examples of such a 
"fine tuning" .On the other hand in the first approximation we neglect an interaction between 
two mesons in the continuum, like DK, etc. 

In present paper we concentrate on the P-wave B, B s mesons and the effects of the 
channel couplin g. W hile the IP levels of the P, D s mesons are now established with good 

or the B,B S mesons only relatively narrow 2 + , 1 + states have been 



accuracy 

recently observed 



271 ]. 281 ] . According to these data the splitting between the 2 + and 1 + 
levels is small, ~ 20 — 10 MeV, while the mass difference between B s (2 + ) and B(2 + ) states 
is again ~ 100 MeV, as for the D s (2 + ) and P(2 + ) mesons. 

The actual position of the B(1P), B S (1P) levels is important for several reasons. Firstly, 
since dynamics of (qb) mesons is very similar to that of qc, the observation of predicted large 
mass shifts of the P s (0 + , 1 + ') levels woul d g ive a strong argument in favour of the decay 
channel mechanism suggested here and in 18j. Secondly, observation of all P-wave states for 
the B,B S mesons could clarify many unclear features of spin-orbit and tensor interactions in 
mesons. Understanding of the decay channel coupling (DCC) mass shifts could become an 
important step in constructing chiral theory of strong decays with emission of one or several 
NG particles. 

The paper is organized as follows. In the next Section we discuss the formalism from 18j , 
extending that to the case of the B, and B s mesons and also to the 1 + states, and discuss 
the mixing between the 1 + and 1 + ' states. In Section 3 the masses of HL mesons, calculated 
in closed-channel approximation, are given. The Section 4 is devoted to the mechanism 
of chiral transitions while in Section 5 our calculations of the mass shifts due to DCC are 
presented. The predictions of the B(J P ), B S (J P ) masses and discussion of our results are 
given in Section 6, while Section 7 contains the Conclusions. In Appendix A a connection 
between the lowest order of ALfcm an d the effective Lagrangian is illustrated. In Appendix 
B the details of our calculations of the masses are given, while in Appendix C the connection 
between FS splittings and the mixing matrix (and angle) of the 1 + states is discussed. 
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II. MIXING OF THE 1+ AND 1+' STATES 



It is well known that in single-channel approximation, due to spin-orbit and tensor interac- 
tions the P-wave multiplet of a HL meson is splitted into four levels with J p = + , lj, 1^, 2 + 



29| | . Here for the 1 + states we use the notation H(L) for the higher(lower) eigenstate of 
the mixing matrix because apriori one cannot say which of them mostly consists of the 
light quark j = 1/2 contribution (see Appendix C). For a HL meson, strongly coupled to a 
nearby decay channel (DC), some member (s) of the P-wave multiplet can be shifted down 
while another not. Just such situationis takes place for the D S (1P) multiplet. 

The scheme of classification, more adopted to a HL meson, in the first approximation 
treats the heavy quark as a static one and therefore the Dirac equation can be used to define 
the light quark levels and wave functions 1Q( |. Starting with the Dirac's P-wave levels, one 
has the states with j = 1/2 and j = 3/2. Since the light quark momentum j and the 
quantum number x are conserved, 1 they run along the following possible values: 
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(2) 



The HL meson w.f. can be expressed in terms of the light quark w.f. - the Dirac bispinors 
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where ^Mrj 2 M 2 are ^ ne corresponding Clebsch-Gordan coefficients. 

Later in the w.f. we neglect possible (very small) mixing between D(ly 2 ),D(l^, 2 ) states 
and also between D s (2^ 2 ), D s (2^ 2 ) states. However, physical P s (l + ) states can be mixed 



1 we use here the standard notation x = =p|j + =| for j 
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via open channels and tensor interaction, while the + and 2 + levels are obtained solely from 
j — | and j — |, respectively. 

The eigenstates, defining the higher 1^ and lower 1^ levels, can be parametrized by 
introducing the mixing angle 0: 

1 3 

|1+) = cos0|j = -) +sin0|j = -), (5) 

and 

1 3 

|1+) = -sin 0| j = -) +cos0|j = -), (6) 

where the mixing angle is defined by the unitary mixing matrix O m i X . In the heavy-quark 
limit the states with j ' = | and j = \ are not mixed, but for finite m,Q they can be mixed and 
definition of the mixing matrix in this basis is rather complicate procedure 10(], which is also 
model-dependent. Therefore it is more convenient to connect the angle <p in ©,([6]) with the 
known factors in the LS basis, where O m i X is well defined in closed-channel approximation 
and factually depends only on the ratio a/t, where a is the spin-orbit and t is the tensor 
splitting. For our analysis we do not need to know details of spin-orbit interaction (see 
Appendix C). 

Then the splittings of the 2 + and + levels are 

M(2+) -M cog = a-0.1t, (7) 

M(0 + ) - M cog = -2a - t, 

while 1^ and 1^- in ([5]) and ([6]) can be expressed through the mixing angle 8 in the expansion 
of these states in the LS basis, where they represent the decomposition of the 3 P\ and l P\ 
states: 

|1+ = cos#| 3 Pi > -sin0| 1 P 1 >, 

|1+ = sin#| 3 Pi > +cos#| l Px > . (8) 

The states 3 Pi and 1 P\ in ([8]) can be expressed through the basis with the eigenstates \ j = |) 
and j = \) |29|: 



■*>-^H> + vfH>' (9) 



-\ —\- = h 

T V3 U 2' 
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TABLE I: The mixing coefficients of the 1 + states for different ratios a/t. 
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Then in closed-channel approximation the following relation can be established between 
angles and 9: 

= -# + 35.264°. (10) 

From ffTUl) it follows that 

1. If 1^ is pure \j — | > state (0 = 0°), then this state is the admixture of the 3 Pi and 
1 P l states with 9 = 35.264°. 

2. If 1^ is pure \j — | > state (0 = 90°), then in the LS basis this state is admixture of 
the 3 Pi and 1 P l states with 9 = -54.736°. 

3. The special case with = —9.74° corresponds to "equal" mixing between the 3 P\ and 
l P\ states with the angle 9 = —45°. 

Therefore the solutions with small |0| < 6° correspond to the mixing angle 9 from the 
range: 29° < 9 < 41° for which the ratio a/t appears to be slightly smaller (larger) unity 
for small negative (positive) 0. Here we consider small negative 0. 

The dependence of the mixing angle on the ratio a/t is illustrated by Tabled taking three 
different ratios a/t. 

Thus 1^ state as the pure \j = \ > corresponds to a = t, while for slightly smaller ratio, 
a/t = 0.95, the admixture of the \j — | > state is ~ 15% and for a/t = 0.917 the admixture 
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TABLE II: The B meson masses (in MeV) (without decay channel coupling) 
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is already 26%. Notice that the physical condition a < t contradicts the heavy-quark limit 
when t — > while a ^ and can have large magnitude. 

The structure of the mixing is important because it defines the order of levels and the 
value of mass shift for the 1 + ' state, as well as the mass shift and the width of another 1 + 
level. It is important that if the coupling to nearby continuum channel is taken into account, 
then as follows from experiment, the position of the 2 + and 1 + levels does not change (within 
1-3 MeV) and just their mass difference A = M(2 + ) — M (1 + ) can be used to define tensor 
splitting: it can be derived that A = (1.25 ± 0.15)i for any a/t. 



III. THE MASSES OF HEAVY-LIGHT MESONS 



In closed-channel approximation the masses of HL mesons, or initial positions of the levels 
without channel coupling), can be calculated in different schemes, e.g. in the LS coupling 



18j, or as in the Dirac type coupling 10j. In Tables [TT1 IHII we give these unperturbed 



masses for the B and B s mesons which are calculated with the use of the relativistic string 
Hamiltonian [f], . In this approach the P-wave masses of HL mesons appear to be smaller 
that in other potential models because they contain negative string corections (see Appendix 
B). 
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TABLE III: The B s meson masses (in MeV) (without decay channel coupling) 
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Calculated masses of the states with j = 3/2, l£ and 2 + , for the ratio a/t = 0.95 (see 
Table [II]) appear to be in good agreement with recent DO Collab. measurements of the B 



meson masses 



281 ] . Such agreement can also be reached for other values a/t close to unity: 



R = - = 1.0 ±0.05, (11) 

and we take the same ratio for the B, B s mesons, and also for the D,D S mesons: for such 
choice the contribution from the j = \ state dominates in the 1^ meson. In particular for 
a/t = 0.95: 

|l+) = 0.9884|i>- 0.1517| ~>, 

and 

|1J> = 0.1517| |> H- 0.9884| |>. (12) 

The masses given in Tables HI] II III are obtained taking the tensor splitting t = 12.2 MeV 
and t = 10 MeV for the B and B s mesons, respectively. The tensor splittings have been 
determined to fit the mass difference M(2 + ) — M(l + ), which has the important property - 
it does not change (within 2 MeV) if DCC is taken into account. 

In Table HV] the masses of the B,B S mesons (from Tables HTllHip are compared to those 
obtained in other models; there the conventional notations 1 + and 1 + ' for the j 1 = | and 



TABLE IV: Theoretical predictions for the B(1P) and B S (1P) masses (in MeV) (without decay 
channel coupling) 
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a Experimcntal data of the DO Collaboration [27] 
^Experimental data from [28] 

j — | states are used. 

Comparison of the masses, given in Table HV] shows that in different papers M(B), M(B S ) 
differ not much, within ±50 MeV, however, the order of levels inside the IP multiplet appears 
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TABLE V: The masses of the D(1P) and D S (1P) mesons in MeV) (without decay channel coupling 
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to be different. In particular, in [J], [jl]] the 2 + level has smaller mass than the 1 + ' while in 
our calculations the 2 + state has always maximal mass. It means that FS parameters a, t, and 
their ratio, as well as the mixing angle (f) between the 1 + and 1 + ' states, can differ essentially 
in given papers. Meanwhile existing experimental limit on the width of _D s i(2536) puts strict 
restrictions on admittable value of the mixing angle <f> (see next Section). 

Finally in Table IVl we give also unperturbed masses of the D(1P) and D S (1P) mesons, 
taking the splitting t = 29 MeV from the mass difference, A = M(2 + ) - M(l+) = 1.31* = 38 
MeV and a/t = 0.95, both for the D and D s mesons. Notice that the position of the P-wave 
mesons does not practically change if a/t = 1.0 (or = 0). 

Given in Table |V] masses show that in closed-channel approximation we have reached 
good agreement with experiment for all D(1P) mesons and also for narrow mesons 
D s (2535), £> s (2573). 

We do not need here to know the details of spin-orbit interaction which at present is not 
ully understood, probably, because of important role of one-loop (or even higher) corrections 



3l| and possible suppression of NP spin-orbit potential observed on the lattice 32J. Here 
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we would like only to notice that in heavy quarkonia the ratios a/t are also close or equal 
unity: 



a/t = 1.04 ±0.08 forx&(LP); a/t = (1.02 ± 0.14) for x 6 (2P), 

a/t = (0.86 ±0.02) for x c (lP). (13) 



IV. CHIRAL TRANSITIONS 



To obtain the mass shift due to DCC effect we use here the chiral Lagrangian (1), which 
includes both effects of confinement (embodied in the string tension) and Chiral Symmetry 
Breaking (CSB) (in Euclidean notations): 

L fc m = i J d 4 x?P + (d + m + M)ip (14) 

with the mass operator M given as a product of the scalar function W(r) and the SU(3) 
flavor octet, 



M = W(r)exp(z l5 ^) : 



where 



(15) 



(16) 



Taking the meson emission to the lowest order, one obtains the quark-pion Lagrangian in 
the form 
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(17) 



Writing the equation (fT71) as AL FC m = — f Vifdt, one obtains the operator matrix ele- 
ment for the transition from the light quark state i (i.e. the initial state i of a HL meson) to 
the continuum state / with the emission of a NG meson ((p a X a ). Thus we are now able to 
write the coupled channel equations, connecting any state of a HL meson to a decay channel 
which contains another HL meson plus a NG meson. 

In the case, when interaction in each channel and also in the transition operator is time- 
independent, one can write following system of equations (see 33J for a review) 

[(H i -E)5 il + V il }G lf = l. (18) 
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Such two-channel system of the equations can be reduced to one equation with additional 



DCC potential, or the Feshbach potential 



(#! - E)G U - V 12 —^—-V 21 G 11 = 1. (19) 
n 2 — & 

Considering a complete set of the states |/) in the decay channel 2 and the set of unperturbed 
states \i) in channel 1, one arrives at the nonlinear equation for the shifted mass E, 

E = E? - ^\ Vn \f)-^—{f\V 2l \i). (20) 
/ E 2 — E 

Here the unperturbed values of E^ are assumed to be known (see Tables [Til EH Ej) , while 
the interaction Uif is defined in (fT7|) . A solution of the nonlinear equation (|20|) yields (in 
general a complex number E = E — ^) one or more roots on all Riemann sheets of the 
complex mass plane. 



V. CALCULATION OF THE DCC SHIFTS 



To calculate explicitly the mass shifts, we will use the Eq. (120]) in the following form: 



<i\V\f> 
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^ £y — m[i] 

where m^°^[i] is the initial mass, m[i] - is the final one, Ef = ujd + ujk is the energy of 
the final state, and the operator V provides the transitions between the channels (see the 
comment after Eq. (JT7j)). 

In our approximation we do not take into account the final state interaction in the DK 
system and neglect the D-meson motion, so the w.f. of the i, /-states are: 

|/ >= V K (p) (g) V D (M f ), \i >= *d s (M,), (22) 

where 

*' f(pH 7& (23) 

is the plane wave describing the A'-meson and ^fr>(Mf), *f> D s (M i ) are the HL meson w.f. at 
rest with the spin projections Mf, Mi, respectively. 
We introduce the following notations: 



wjr(p) = \/p 2 + m 2 K , u D (p) = \/p 2 + m 2 D , (24) 
13 



so that in the final state the total energy is Ef = ujd + uk, while 

Tf = Ef — mo — rriK (25) 

is the kinetic energy. Also it is convenient to define other masses with respect to nearby 
threshold: m t hr — m K + m D, 

E = m (0) [D s ] -m D - m K , 5m = m[D s ] - m (0) [D s ] , (26) 

A = Eq + 5m = m[D s ] — mo — mx, (27) 

where A determines the deviation of the D s meson mass from the threshold. In what follows 
we consider unperturbed masses mo(J p ) of the (Qq) levels as given (our results do not change 
if we slightly vary their position, in this way the analysis is actually model-independent). 
Using these notations, the Eq.(l20l) can be rewritten as 

Eo-A = F(A), (28) 



where 

def 



2 



^ (Mi V p, M f 



and 



Mi 



V" 



e 



p,MA = - / ( Mi )a\r\ l5 ?^* D (M f ) ^==d\ (30) 
1 J Jk ^2u K (p) 

The function ^"(A) for negative A diminishes monotonously so there exists a final (crit- 
ical) point, 

£ crit = JF(-O). (31) 

Thus, while solving the Eq.( l28l) . one has two possible situations: Eq < E^ A and Eq > Eq* 11 . 

In the first case Eq. fTSBl) has a negative real root A < (see Fig. [T|) and the resulting 
mass of the D s meson appears to be under the threshold. In the second case Eq. (!28|) has 
a complex root A = A' + iA" with positive real part A' > (see Fig. [2]) and negative 
imaginary part A" < 0. To find latter solutions one should make analytic continuation of 
the solution(s) from the upper halfplane of A under the cut, which starts at the threshold, 
to the lower halfplane (second sheet). This solution can be also obtained by deforming the 
integration contour in Tf(p). In actual calculations we take infinitesimal imaginary part A", 
proving that A does not change much for finite A" (the similar procedure has been used in 
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FIG. 1: Eq.([2S} for the case E < Eff 



^ E -A 



-100 -80 - 60 -40 -20 

A, MeV 



FIG. 2: Eq.(T28]) for the case E > Eg 
'\Eo-A 



f (A) 



20 3C 
A, MeV 



[l^D- Finally, the resulting mass of the _D S meson proves to be in the complex plane at the 
position A' — i|A"|, i.e. the meson has the finite width T = 2A". 

For further calculations we should insert the explicit meson w.f. to the matrix element 
(!30|) . As discussed above, in a HL meson we consider a light quark q moving in the static 
field of a heavy antiquark Q, and therefore its w.f. can be taken as the Dirac bispinor: 



g(r)ttjiM 

(-i)^/(r)n^ 



I {f + g 2 ) r'dr = 1, 



(32) 



where the functions g(r) and f(r) are the solutions of the Dirac equation: 



1 + x 



9 



g-(E q + m q + U-V c )f = 0, 



r 

f + ^—^f + (E q -m q -U -V c )g = 0. 



(33) 
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Here the interaction between the quark and the antiquark is described by a sum 
scalar potential and the vector Coulomb potential with a s = const: 



-, p = -a„. 
r 6 



U = ar, V c = 
Introducing new dimensionless variables 

x = ta/o 7 , e q = E q / y/a, \i q 
and new dimensionless functions 



m, 



o. 



g = a 



3/4 



G(x) 



x 



f = <r 



3/4 



F(x) 



x 



OO 

J {F 2 + G 2 ) dx = l, 



we come to the following system of equations: 



G' + —G - f e q + n q + x + - J F = 0, 

X \ X J 



F>-*F 

X 







~ q — Hg — x H — J G = 0. 



This system has been solved numerically. 
Using the parameters from the papers 



a = 0.18 GeV 2 , a, = 0.39, 



210 MeV, m q = 4 MeV, 



we obtain the following Dirac eigenvalues e: 





Qq, Hg = 0.01 


Qs, y, a = 0.5 


-1 


1.0026 


1.28944 


+1 


1.7829 


2.08607 


-2 


1.7545 


2.08475 



and corresponding eigenf unctions G, F are given in Figs. [HI HJ) . 
Later we use the simplified notations for the quark bispinors 

dcf ,i,0,M 2 



def , i,l,Mi 
1 = V* ' 



I I A/r \ dcf / f.i.-M'a 

^ 3 (M 3 ) = VI 



Now, using explicit expressions for the spherical spinors 

/ j+M Y 

^{+1/2,I,M = r > ^-1/ 



/2 



2,LM 



j-M+1 v 



/2 



i+M+i v 

2j+2 r ',Af+l/ 
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FIG. 3: Gi i 2,3(x) functions 




FIG. 4: Fi )2 ,3(a;) functions 
i(x) 

:F 2 (x) 




and the expansion : 



l.M ^ P ' 



(42) 



after cumbersome transformations (which are omitted in the text) we obtain the transition 
matrix elements: 



Vl2 









V2 



D ipr 



4(M,) a|r| 75 — ^ 2 (M 2 ) d 3 r 

Ik v /2u; K (p) 



fKy/VK{P) \V cr 



\ ~J ' ( 43 ) 
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FIG. 5: $0,2(5) functions 




*o(q) 



v 32 


M 3 ,M2 j 







3 ipr 



^(M 3 ) ( r|r|75^V2(M 2 )— J== 
Jk y/2u K {p) 



d 3 r 



^ $ ( P \ J 47T Y> 



where 



-1 +2 
-2 +1 



*o(g) = / jo{qx)xdx Gi{x)F 2 (x) - F 1 (x)G 2 {x) 



^2(9) = / 32{qx)xdx G 3 (x)F 2 (x) - F 3 (x)G 2 (x) 



(44) 



(45) 



Notice that because of different signs of the Fi(x) and F 2j 3(x) functions (while the £1,2,3 
functions are all positive) on almost all real axis, the integral $ 2 appears to be strongly sup- 
pressed in comparison with the integral $o- This fact is confirmed by numerical simulations 
(see Fig. EJ). 

Finally, introducing universal functions 



^o, 2 (A) 



$2 (pSj!) 

a f p(T f )uj D (T f )dT f °H vW 



2^f 2 K J T f + m D + m K T f - A 


r (T\ M T /WT/) 2 f p( T f) \ 

J- 0,2 U f J = — 73- ' 77T— : ' ^0,2 — 1=- 

71 Jk T f + m D + m K ' \ y/a J 



(46) 
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we come to the following equations to determine meson masses and widths: 



D s (0+) 
D.(lt) 



D s {2 



3/2/ 



E [0 + ]- A = ^o(A), 

Eb[l£] - A = cos 2 ■ ^o(A) + sin 2 • F 2 (A), 

£b[lj] " A' = sin 2 • ^o(A') + cos 2 • JS(A') 

r[l+] = sin 2 • fo(A') + cos 2 • f 2 (A'), 
3 



E [2 



, - A' 

3/2J " 5 



r[2+ 2 ] = -.f a (A'). 



(47) 



VI. RESULTS AND DISCUSSION 

In this chapter, using the expressions (14 7p to define the D s and B s meson mass shifts, we 
present and discuss our results. We will take into account the following pairs of mesons in 
coupled channels (i refers to first (initial) channel, while / refers to second (decay) one): 



i 


D s (0+) 


Ds(l + ) 




f 


D(0-)+K(0~) 


D*(l-) + K(0-) 


D*(l-) + K(0-) 


i 


B s (0 + ) 


B s (l + ) 


B s (2+) 


f 


B{Q-) + K{0~) 


B*(l-)+K(Q-) 


B*(l-) + K(Q-) 



(48) 



In our calculations we use here the following meson masses and thresholds (in MeV): 



m D + = 1869, m D + + m K - = 2363, 

rriD*+ = 2010, mo*+ + Trig- = 2504, 

m B + = 5279, m B + + m K - = 5772, 

m B * = 5325, rriB* + m K - = 5819. 



(49) 



Villi A priori one cannot say 



The results of our calculation are presented in Tables IV 
whether the \j = |) and \j = |) states are mixed or not. In 
mixing at all, the width T(D s i(2536)) = 0.3 MeV is obtained, while now the experimental 



35j in the case when there is no 



limit is T < 2.3 MeV 26] and recently in , 



the width T = 1.0 ± 0.17 MeV has been 



measured. Therefore small mixing is not excluded and here we take the mixing angle 
slightly deviated from = 0° (when there is no mixing at all). Then one can define those 
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TABLE VI: D s (0 + )-meson mass shift due to the DK decay channel and S s (0 + )-meson mass shift 
due to the BK decay channel (all in MeV) 



state 




m (theor) 


m (cxp) 


5m 


D s (0+) 


2467 


2331 


2317 


-136 


B s (0+) 


5805 


5700 


not seen 


-105 



TABLE VII: The D s (l + ), D s (2 + ) meson mass shifts and widths due to the D*K decay channel 
for the mixing angle 4° (all in MeV) 



state 


m (o) 


m (theor) 


m (exp) 


-p(theor) 
1 (D*K) 


-p(exp) 
1 (D*K) 


5m 




2550 


2440 


2460 


X 


X 


-110 




2537 


2535 


2535 


1.1 


< 1.3 


-2 


D s (2+ /2 ) 


2575 


2573 


2573 


0.03 


not seen 


-2 



angles (f> which are compatible with experimental data for the masses and widths of both 1 + 
states. The limiting angle \(f>\ = 5.7° (given in Table HX|) corresponds to the mixing between 
the 3 Pi and x Pi with 9 = 41° in LS scheme. 

The value £ = (0.995) 2 for the 1^ states provides large mass shift (~ 100 MeV) of this 
level and at the same time does not produce the mass shift of the l£ level, which is almost 
pure j ' = | state. For illustration we show the scheme of the 1 + , 2 + shifts on Figures [61 [71 
We would like to stress here that the dependence of the shift on the heavy non-strange meson 

TABLE VIII: The B s (l + ), B s (2 + ) meson mass shifts and widths due to the B*K decay channel 
for the mixing angle 4° (all in MeV) 



state 


m (0) 


m (theor) 


m (exp) 


-p(theor) 
1 (B*K) 


p(exp) 
(B*K) 


5m 




5835 


5727 


not seen 


X 


X 


-108 


B.ot) 


5830 


5828 


5829 


0.8 


< 2.3 


-2 


B s {^l/ 2 ) 


5842 


5840 


5840 


< 10" 3 


not seen 


-2 
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TABLE IX: The mixing coefficients in Eq.flU} for the DCC shifts of the B s , D s mesons = 5.7°). 









1J 




a/t 


cos 2 (j) 


sin 2 


D S (1P),B S (1P) 


0.96 


(0.995) 2 


(0.100) 2 


D S (1P),B S (1P) 


1.0 


1.0 






FIG. 6: Scheme of D S {1 + ,2 + ) shifts due to chiral coupling 




2600 



2550 



2500 



2450 



2400 



mass (or, equivalently, on the heavy quark mass) is rather weak, which follows directly from 
the Eq. (1471) via inverse mass expansion. 

If one uses one-gluon exchange potential for tensor interaction, then to obtain the splitting 
t ~ 12 MeV for the B and B s mesons and t = 29 MeV for the D and D s mesons one needs 
to take a s ([iFs) ~ 0.39 for all P-wave HL mesons. However, in OGE approximation the 
spin-orbit splitting does not satisfy the condition a = 0.95£ for this value of a s . This 
fact possibly indicates on important role of one-loop or even higher radiative corrections, 
observed in heavy quarkonia 3j, and also on p ossible ~ 20% suppression of NP spin-orbit 
confining potential observed on the lattice 32 ] . 
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FIG. 7: Scheme of l? s (l + ,2 + ) shifts due to chiral coupling 



B.(1S) 



B,(l£) 



B s (23/2) 



5900 



5850 



5800 



5750 



5700 




VII. CONCLUSIONS 

We have studied the mass shifts of the D a (0 + , 1+ ) and B s (0 + , 1+ ) mesons due to strong 
coupling to the decay channels DK, D*K and BK, B*K. To this end the chiral quark-pion 
Lagrangian without fitting parameters has been used. 

We have shown that the emission of a NG meson, accompanied with the 75 factor, gives 
rise to maximal overlapping between the higher component with j = | of the P-wave meson 
(D s , B s ) bispinor w.f. and the lower component (also with j = |) of the S'-wave meson w.f. 
in considered S'-wave decay channel. Due to this effect, while taking the w.f. of the IP 
and IS* states with the use of the Dirac equation, large mass shifts of the + , 1 + ' states are 
obtained. 

The widths of _D s i(2536) and _B s i(5830) are also calculated. To satisfy the experimental 
condition r(D s i(2536)) < 2.3 MeV the following limit on the mixing angle (between the 
\j — I > and \j — ~ > states) is obtained, \(f>\ < 6°. This restriction implies that the mixing 
angle 9 between the 3 Pi and l P\ states in LS basis lies in the range: 29° < 9 < 41°. For 
given angles 9 the ratio of the spin-orbit and tensor splittings appears to be close to unity 
a/t = 1.0 ± 0.05, as it takes place for the Xb an d Xc mesons. The value of tensor splitting 
t has been defined from the mass difference, M(2 + ) — M(l + ), which is not affected by the 
coupling to the decay channel. 

Calculated masses of the 2 + and 1 + states are in good agreement with experiment for all 
D, D s , B, B s mesons. 

For the + , 1 + states, taking into account the mass shifts, the following masses are 



22 



predicted: 

• M(B*) = 5695(10) MeV which coincides with M(B*) = 5695(10) MeV, 

• M(B S (1+')) = 5730(10) MeV, close to M(B 1 (1+')) = 5732(10) MeV. 
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APPENDIX A: CONNECTION BETWEEN THE CHIRAL QUARK-PION LA- 
GRANGIAN AND THE EFFECTIVE CHIRAL LAGRANGIAN 



The interaction of pions with quarks was introduced and developed in 21], 22J, see 17] 
for recent applications. The effective chiral Lagrangian AL e // contains one new parameter 
9% 

L« q = |f (Al) 

and has the form of pseudovector coupling, known from phenomenological applications in 
the pion-nucleon systems. As was argued in [2lj], g q A at large N c tends to unity. 

In Eq.(17) we have used another form of the quark-pion interaction, derived directly from 
the QCD Lagrangian in [20I and not containing new parameters, 

A4cAf = / *P(x)a\x\ l5 ^^(x)d 4 x. (A2) 

In {^(J the connection between (1A1|) and (1A2j) was established and here we repeat the deriva- 
tion for the convenience of readers. 

Consider application of (1A2I) to the case of pionic transition between states 4>i(x) and 
ip2{%) of the quark in heavy- light meson. Dirac equations for ipi(x) can be written as 

(ctp + f3(m + a\x\) + V cou i)ip 1 = (A3) 

^(-«P + (3(m + <t|x|) + Vcoui) = e 2 ip 2 - (A4) 
Expressing in ( 1A3I) . ( 1A4I) the term cr|x|7 5 via ap, [3m etc. and summing two equations, one 



gets 



al fcm = ^p^2(-2m7 5 7T + - £i)tt + 7 B /3ap7r)^i. (A5) 

^ J IT 
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Since ji = —i{3cti, and (£2 — £i)tt = i§t^{t) 7r(t) ~ e i( - e2 one can rewrite the last two 
terms in (1A5I) as 7^75^71", and finally one arrives at 

Al fcm = ^V^2(-2m7 5 7r + 7 M 7 5 d M 7r)V>i. (A6) 

Comparing ( 1A1I) and (1A6I) . one can see that in the chiral limit, m q — » 0, two expressions 
coincide. However, for nonzero m, e.g. for strange quark having the mass m s ~ 0.2 GeV at 
low scale ~ 1 GeV first term in ( 1A6I) is becoming essential. Moreover, our expression 
( 1A2[) is only the first term in the expansion of the exponent (15) in powers of the pion field, 
and therefore this general Lagrangian can be used for decay channels with the production 
of two or several pions. 



APPENDIX B: MASSES OF HEAVY-LIGHT MESONS 



To calculate masses and different matrix elements (m.e.) of a HL meson (qb, qc, or qb) 
we use here the relativistic string Hamiltonian H u , derived in 24|. For this Hamiltonian the 
spin-averaged mass M cog (nL) is given by simple formula: 



M cog {nL) = M (nL) + A SE + A 



str j 



(Bl) 



where M is the eigenvalue (e.v.) of the spin-independent part H of the Hamiltonian H w , 
which coincides with well-known spinless Salpeter Hamiltonian (SSH): 



p 2 + m 2 b + V (r), (B2) 

H^ndr) = M oVnl (r). (B3) 

However, the mass (IBlj) contains negative (string) correction and therefore in our approach 
for a given static potential the levels with L 7^ lie lower than for SSH. Also the mass (1B1I) 
does not contain an overall fitting constant but takes into account NP self-energy term A^^ 
for a light quark ( which is calculated explicitly in 38j|). 
The static potential 



Vn(r) 



ar 



4a B (r) 



(B4) 



is taken here from 



251 ] with the vector oupling a_e(r) for rif = 3. The solutions of flB3j) 



define M (nL) and m.e., in particular, 



u q (nL) = (yV + m? q ) nL , u Q (nL) = (yV + m 2 Q ) nL , 



(B5) 
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TABLE X: The constituent masses io q and cj& (in MeV) for the B{\P) and B S (1P) mesons (m u ^ = 
0, m s = 200 MeV, m b = 4780 MeV) 





B(1P) meson 


B S (1P) meson 


u q (lP) 


680 


730 


u b (lP) 


4836 


4840 


^red 


598 


634 



which apear to be the dynamical (constituent) mass of a light quark uo q and ujq for a heavy 
quark. Their values for the B and B s mesons are given in Table 1X1 together with the reduced 

UlqUJb 



mass: oVed 



U)q+U) b ' 



As seen from Table IXl the kinetic energy of a light (strange) quark u> q (lP) are not small 
and this fact is important for the fine structure analysis. 

In the mass formula (1B1I) the correction A$e comes from NP self-energy contribution 
(which is equal zero for the b quark), taken here in the simplest form when self-energy 



contribution of the c quark (< —15 MeV) can be neglected, because it is sma 



1 as compared 



to the pole c-quark mass known at present with the accuracy ±100 MeV 26|). For a light 
quark A S e has been defined in {38]: 

1.5(777, 



*SE 



nL) 



(B6) 



in which small correction, ~ 3 4- 6 MeV (defined in 25J) is neglected. The factor r] u ^) = 1.0 
for a light quark and rj s = 0.65 for the s quark with m s = 220 MeV. 

The string correction A str for the IP -wave B(B S ) mesons is equal A str ~ — 27(— 21) 
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TABLE XI: The masses M ,M cog (lP) and A SS (1P), A str (lP) (in GeV) for the B,B S mesons 
(m s = 220 MeV, m u(d) = 0, m b = 4780 MeV) 







B.(1P) 


M (IP) 


5885 


5925 


A 5 £ 


-126 


-70 


A s j r 


-27 


-20 


M C09 (1P) 


5.732 


5.835 



MeV 



25] . This negative contribution to M cog improves an agreement with the experimental 



masses of B(2+) and 5(1+) mesons |28|. In TableEDthe eigenvalues M (1P) and M cog (lP) 
together with As^(lP) and A str (lP) are given. 

The difference between the e.v. M (B) and M (B S ) is only ~ 50 MeV for m s = 200MeV, 
so that additional 50 MeV difference in M cog (lP) for the B and B s comes from the self-energy 
terms. 



APPENDIX C: FINE STRUCTURE SPLITTINGS 



To define FS splittings of a HL meson we follow here the approach, where the LS basis 



is used and the analysis of FS can be done in general terms 



. If one introduces tensor 



splittings t(nP) and spin-orbit splitting a{nP) then the masses of the 2 + state (j ; = |) and 
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+ state (J = |) can be written as 

M(2 + ) = M cog + a - O.lt, (CI) 

M(0 + ) = M cog -2a-t, (C2) 

while the 1 + states, 3 Pi and 1 Pi, are mixed. 

The mixing matrix can be expressed through the splittings a and t: 

a -l t -^t 

Orrn* = | J ' _ 6 | . (C3) 

Then the eigenvalues and eigenvectors of this matrix define "higher" and "lower" masses 
Mn and Ml with J p = 1 + and the decomposition of their w.f.. The mass splittings are 

Mh = M cog - I(2a - t) + ^ v / (2a-t) 2 + 32(a-t) 2 ; (C4) 

M L = M cog - i(2a - t) - ^ v / (2a-t) 2 + 32(a-t) 2 . (C5) 

Each of these levels is a decomposition of the Z P\ and : Pi states. From (1C3|) it is evident 
that the weights in those decompositions depend only on the ratio 

R=- t . (C6) 

Just the value of this ratio defines the order of levels inside the nP multiplet and for given 
a/t the mixing angle 9 in (9) can be easily calculated. With the use of the relation (10) the 
connection between the mixing angle <f> in the basis \j = |), \j = |) (see the definition in 
(6), (7)) and the mixing angle 9 in the LS scheme, for which 

|1+) = cosfl| 3 Pi) -sin0| 1 P 1 ), 

have been established: <ft = —9 + 35.264°. 

The value a/t = 1.0±0.05 provides large mass shifts of the D s (ljj) and P s (l^-) levels and 
at the same time keeps the position of the D S {1\) and B s (l^) unchanged (with accuracy 
2 MeV). Notice, that for all multiplets with a/t = 0.95 the mass difference between two 
narrow levels is 

M(2 + ) - M{±1) = 1.31* for a/t = 0.95, (C7) 
M(2 + ) - M(l+) = 1.40t for a/t = 1.0, 
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so that this relation can be used to define the parameter t from experiment. 

To interpret this splitting t one can use well-known perturbative expression, taking one- 
gluon-exchange interaction and neglecting higher in a s corrections, which however may be 



important 



3l|: 



t{nP) = j^^K (C8) 

3 LO q UjQ 

Then for cxfs — 0.39 one obtains the values: t = 11.7 MeV for the B s mesons ((r~ 3 )ip = 
0.080 GeV 3 ) and t = 12.1 for the B(1P) mesons ((r~ 3 ) 1P = 0.0765, which are close to 
those, used in our analysis. For the D(1P) and D S (1P) meson the same otps — 0.39 gives 
t D = t Ds = 29 MeV ((r- 3 ) D = 0.052 and (r~ 3 ) Da = 0.055 GeV 3 ). 

However, in OGE approximation the situation with SO splitting is not so simple. This 



splitting can be presented in the convenient form from 
proportional tHq 2 : 



29| (here we also keep the term 



where the factor 



<nP) = ( JL + _L | A (nP) + t(nP), (C9) 



A(nP) = ^a FS (r- 3 ) nP - a^- 1 )^. (C10) 



To satisfy the condition a = 0.95t the term A(1P) has to be small as compared with t. 
However, in OGE approximation and for linear confinement this term in (1C10I) appears to 
be negative and not small for olfs ~ 0.39. The reason for that needs a _special investigation 



and can mean that either higher radiative corrections are important 31], or a suppression 



of confining potential in spin-orbit term ,observed on the lattice [32| , is essential. 

Notice that the Coulomb-type order of levels, i.e. M(0+) < M(l+) < M(l+) < M(2+), 
takes place only for not small ratio | > 0.606. In our case with | = 0.95 this condition is 



satisfied and the level 1^ lies below the 2 + level, 



M H (l + ) < M(2+). (Cll) 

Just this order of levels is observed in the D(1P) multiplet where the central value of the 

261 ] . The FS splittings of the P-wave HL 



wide 1 + ' level is smaller than the mass of 2 + state 



mesons with a/t = 0.95 are given in Tables HTj IIIIl IVl 
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